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Abstract 

The interaction of light quarks and instanton liquid is analyzed at finite density of 
quark/baryon matter and in the phase of nonzero values of diquark (colour) condensate. It 
is shown that instanton liquid perturbation produced by such an interaction results in an essen- 
' tial increase of the critical value of quark chemical potential \i c which provokes the perceptible 

increase of quark matter density around the expected onset of the colour superconductivity 
phase. 
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The intriguing history of studying the structure of strongly interacting matter abounded in many 
! interesting but somewhat contradictory phenomenological results up to the advent of quantum chro- 
modynamics (QCD). Since the time of its discovery the possible existence of quark-gluon plasma 
(QGP) phase is one of the most fundamental predictions of QCD. This phase is phenomenologically 
understood as a matter state where quarks and gluons being under extreme conditions of high tem- 
perature (T) and(or) quark/baryon density (driven by the corresponding chemical potential fi) are 
^ ! able to propagate relatively freely over the considerable (of course, on the hadron scale) distances, 
j Even the simplest quantitative estimates of QGP physical characteristics were so verisimilar that they 
\ allowed to launch the process of QGP experimental study in ultrarelativistic heavy ion collisions. 
5^ ' And although the first series of experiments at SPS CERN and AGS BNL have provided us with 
encouraging results a new generation of experiments is necessary to make convincing conclusions. 
The physical programme for new experimental research certainly needs serious quantitative analysis 
of the QCD phase diagram. 

It seemed lattice QCD could be able to resolve the problem starting even on the first theoretical 
principles. However, up to now this approach was effective and apparently reliable in producing the 
results mainly around the temperature axis of the (/i — T)-plane and, unfortunately, relies on the field 
theory formulation in imaginary time. Nevertheless, lattice QCD has corroborated the phenomeno- 
logical expectations that its phase structure is quite rich |l]. In particular, nowadays we know that 
the chiral symmetry of initial QCD Lagrangian becomes broken and the quarks are getting to be 
confined when the temperature decreases below a critical value. Moreover, the lattice calculations 
show us that both physical phenomena result from the corresponding phase transitions and their 
critical temperatures practically coincide. The latter observation indicates that both fundamental 
phenomena could be tightly knitted and it should heuristically play an important role to illuminate 
the confinement mechanism. 

The features of the QCD matter at high baryon densities remains poorly understood in spite 
of the recent splash of renewed activity [J2J. One of the major theoretical reasons comes from se- 
rious technical difficulties of lattice simulations with Monte-Carlo techniques due to the complex 
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fermion determinant in QCD at finite quark densities. From the phenomenological side the recent 
experimental data are not very informative for developing the theory along this direction. It is clear 
that present and future experiments with heavy ions at higher and higher energies should produce 
the strong interacting matter rather at higher temperature and closer to the temperature axis of 
the (/i — T)-plane. In such a situation the phenomenological treatment of studying along the /i-axis 
should basically rely on the astrophysical observations of the compact stars. In the meantime, it was 
demonstrated that the diquark pairing in colour anti-triplet attractive channel induced by the instan- 
ton interaction vertex should be much more effective than an one-gluon exchange and at sufficiently 
large quark/baryon densities and small T could lead to colour superconductivity and a sophisticated 
picture of the QCD phase diagram f2|. 

The present paper continues our preceding one || where the light quark interactions with the 
instanton liquid (IL) at nonzero values of quark/baryon chemical potentials were investigated in the 
phase of broken chiral symmetry within the improved method of calculating the generating function 
of the IL theory |4| . The new treatment of the functional integral is mainly related to an idea t take 
into account the IL back-reaction upon the quark presence which was always considered negligible. 
In fact, this effect is pretty weak and should manifest itself in the leading order of expansion in the 
effective coupling. Nevertheless, it was qualitatively argued that the interaction of IL and quarks 
could increase the quark matter density around the onset of expected colour superconductivity. As 
known ||, || this density occurred to be surprisingly small (n q ~ 0.062 fm~ 3 ) and, hence, the diquark 
phase would come to play already at the density of nuclear matter (n n ~ 0.45 fm -3 and quark matter 
density is taken to be three times larger than the conventional nuclear one). 

Here we obtain the quantitative estimate of the corresponding critical chemical potential p c ana- 
lyzing the system of the Gorkov equations for the colour superconductor which we further calculate. 
We do not introduce any simplification of the initial Lagrangian trying to deal with the 'exact' four- 
quark interaction which is generated by (anti-)instantons. Such an intention is simply compulsory 
because it is dictated by the precision of problem where the effects of the IL perturbation are expected 
to be of the same order of magnitude compared to the terms which could usually be ignored. 

Let us point out that in the IL approach [|7| the generating functional has the following factorized 
form 



Here the first factor provides us with information on the gluon condensate whereas the fermion factor 
serves to describe the quark practice in the instanton environment ||, [0]. In what follows we 
use the notations of Ref. || where the dimensionless variables (motivated by the form of interquark 
interaction kernel) were introduced, for example, for the chemical potential it looks like p — > pp/2 
and for the momenta kip/2 — > ki, i = 1, . . . ,4 where p is the average size of pseudoparticle (PP). 
The quark determinant Z^ may be transformed with the auxiliary integration over the parameter A 
to the following form (for the colour SU(3) group with the quarks of two flavours Nf = 2): 



and N is the normalizing factor. For clarity we take it to be equal a unity but, in principle, it 
could play a role of free model parameter. This factor is inessential for the models with the fixed 
value of the packing fraction parameter np 4 but in the model where it is admissible for variation the 
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weak logarithmic dependence on N appears. The factors 2 in Eq. ([!]) come with making use of the 
dimensionless variables. The term of four-fermion quark interaction V may be directly expressed 
with the chiral components as 
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it being known that for the right hand fields the substitution L — > R should be performed. The 
kernels Lj?/ are defined by the functions hi, % = 1, . . . ,4 and the zero modes (the solutions of the 
Dirac equation with the chemical potential p in the PP field) 

h 4 (k A , k; p) = jr{(k - P - ik 4 )[(2k 4 + ip)f x + i{k - p - ih)f 2 ] + 
+ (k + p + ik 4 )[(2k 4 + i^i) ft -i(k + p + ih)f+]} , 

hi(k 4 ,k;p) = '^—j j(2& - p)(k — p - ik 4 )fi + (2k + p)(k + p + ik 4 )ft+ 
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2(k - p)(k - p - ik 4 ) - -(p + ik 4 )[kl + (p — p) 2 ] 
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2(k + fi)(k + fi + ik 4 ) + r(p + ifa) \k\ + (p + fi) 2 } 
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where k = \k\ if the spatial components of 4- vector k u are considered and 
f± _ hiz^Koiz*) - Ioiz^K^z*) f ±_h(z ± )K 1 (z ± ) ± _p 



ff= 1V 7 uv ' r K 7 1V ' ,ft= 1V J r ' , z ± = ^ki + (k± f iy 
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with the modified Bessel functions It, Ki (i = 0, 1). Let us introduce the scalar function h(k 4 , k; fi) re- 
f- 
lated to three-dimensional component hi(k 4 , k; fi) = h(k 4 , k; p) , i — 1, 2, 3 (we omit the arguments 

of functions hi when it does not mislead). 

L^(p,q;p) = S ik (p;p)e kl U? U}?e™S+(q; -p) , 
* 

S(p; p) = (p + ip)~ h + (p; p) , S + (p; —p) =h~(p] —p)(p + ip) + where it is valid for the conjugated 

function hp (p\ — p) = h^(p; p) and e is an antisymmetric matrix with e X2 = — e 2i = 1- Here p ± and 
other similar designations are used for the four- vectors spanned by cr^-matrices, = (±z<r, 1), (a 
is the three- vector of the Pauli matrices), for example, p = p v ot an d U is a matrix of rotations in 
the colour space. The similar relations are valid for the right hand components 

(Vf R f Vf ) = C,(p/) *2#(p/,?/;a0 V>f'*fa) > 



with the kernel 



^S(p, t, /*) = r k ( v - p)e kl u? ojv^zjfa -p) 



where T(p; p) = (p + ip) + h~(p; p) , T + (p; —p) =h + (p', —p)(p + ip)~ ■ The components of matrices 
(p + ip)^ 1 and h^(p; p) commute because the vector-function h(p) is spanned by the vector p only. 
Then the following identities are easily understood 

T(p; p) = S + (p; -p) , T + (p; -p) = S(p; p) . 

In what follows we omit the //-dependence of matrices T, T + as the chemical potential enters the 
matrix T being always positive and the matrix T + being negative only. Besides, two other identities 
would also be helpful 

a 2 T T (p)a 2 = T + (p) , a 2 T +T (p)a 2 = T(p) , 
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where T T means a transposed matrix. 

Dealing with the following averages interesting to study the diquark condensates 



(V>i L a? (p)^?(?)) = e i2 71-4 <*(p + ?) > 

(V^ot(p)^S(ff)> = <W * A *(p ~ ?) Gi*{P) , 
and making use of the effective action of Eq. (|I|) one may obtain the Gorkov equations similar to 
ones treated in Ref. [HI 
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F H (p) - H L (p) G RL T {-p) = , 



(2) 



G RL (p) - E h (p) F +L T (p) = 1 



+L T 1 



where 



Goip) 



= —2 (p + ifi)^ means the free Green function, E R (p) = A R T(p) e T T (—p), 
E L (p) = A L T + (p) e T +T (—p) and A L ' R denotes a gap. The form of E-matrices results from the 
kernel structure of equations if averaging over colour orientations done (remember that we consider 
colour stochastic ensemble). In order to complete Eqs. (Q) we need the following gap equations 
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T+(q) F (q) T +T (-q) - T + (-q) F T (q) T (q) 



N C (N C 



dq 
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T(q) F\q) T 1 \—q) - T(-q) F L T (q) T'\q) 



Apparently, the right hand sides of these equations are proportional e because they are the differences 
a matrix and its transposed form. Similar equations are valid for the conjugated matrices 



F +LT {V) [Goip) 

G RL (p) [G (p) 

F +RT {p) [G+(p) 

G LR {p) \G+(p) 
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-p) E +/? (p) = 



-F R (p) E+*(p) = 1 , 

_ Q LR T { _ p) S+ L (p) = ( 

-F\p) E +i (p) = 1 , 



(3) 



with the corresponding gap equations 
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T T (—q) F+ R T (q) T(q) - T T (q) F +R (q) T(-q) 
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T +T (-q) F+ L T (q) T + {q) - T +T (q) F +L (q) T + {-q) 



where H +R (p) = A +R T +T (—p) e T + (p) and E +i (p) = A +L T T (—p) e T(p). In this paper we limit 
ourselves to treating the diquark condensate only. However, as known the mixed phase of nonzero 
values both chiral and colour condensates could exist at \i c ~ 300 MeV realizing the transitional 
regime for the onset of colour superconducting phase. In order to bring this phase to the play the 
equation system should be extended including another averages as 

W/SMltfb)) = Sf 9 <W ^ s(p - q) G^ RR ( P ) . 
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From the Eqs. (0) and fl3|) we find 

G LR ( P ) = G+(p) + G+(p)z R ( P ) f +rt ( p ) , 

F +RT (p) = G LRT (-p) S +L (p) G+(p) . 

With the auxiliary matrices C +L (p) = A +L T + (— p) T(p), C R (p) = A R T(p) T + (—p) we may rewrite 
the matrices £ as 

Z+ L (p) = e C +L (p) , E R (p) = C R (p) e . 
The identities for T-matrix mentioned above help to show the validity of the relations 

e C+ L T {-p) e T = C +L {p) , 

e C R T (-p) e T = C R (p) . 

And for the Green function we have 

G LR (p) = G+(p) + G+(P) C R {p) e G LR T {-p) e C +L {p) G+(p) . 

Combining the properties of C matrices together with the identities for free Green functions 

°2 G± T (p) a 2 = Gf (p) , 

we are able to get the complete equation for calculating the function G LR in the form 

e G LR (-p) e T = Go (-p) + G (-p) C +L {p) G LR {p) C R {p) G (-p) . 

The vector- function h(p) being spanned by the vector p helps to conclude that all the matrices Gq , 
C +L , qr com mute with each other. Then searching the solution for the Green function G LR by 
iterating one finds immediately that the Green function commutes with those matrices and finally 
obtains the equation to calculate it in the form 

G LR (p) = [1 + H{p)\ G+{p) + H\p) G LR {p) , 

or 

[1-H(p)\ G LR (p)=G+( P ) , 

where H(p) = G$(p) C R (p) G (-p) C +L (p). 

The structure of matrices H(p) allows us to conclude that their sum 

H{p) + H(-p) = aip) 

and product 

H{p) H{-p) = ^p) , 

are proportional to the unity matrices (it is clear by definition that at{—p) = aip), fti—p) = flip))- 

Introducing notation g v = h v {—p) we are able to present the functions aip) and flip) in the 
compact form as 

aip) = 4 A* A +L [-4 Aip) (hg) + 2 (p 2 + /?) (^) (/)] , 

where 

Aip) = ip 2 + /j 2 ) (hg) - 2i n p(g±h - Kg) . 

and a scalar product is naturally defined (hg) = J2t=i hidi together with the functions (/i 2 ) and ig 2 ) 
squared. Getting the last term of A(p) we used the designation of scalar function h mentioned above 
and for the function fl we have 

• 2 '- 2\2 /„2\2 



flip) = lQ (A R Y (A+ L Y (p 2 + ^ 2 Y ih 2 ) 2 ig 2 ) 



As both functions a and (3 are spanned by the unit matrices the solution for the Green function 
may be given in 



c lr m = Gt(p) [1 ~ H(-p)} 
[P) - l- a (p)+(3(p) 



The gap equation then looks like 



l 2A rdp a(p)-2/3(p) 



A N C (N C - 1) J 7T 4 A+ L (l - a(p) + (3(p) ' (4) 

We are interested in the solution of the form A R = A +L = A L = A +R at A < 0, which is dictated by 
the symmetries of four-quark interaction kernels. Let us remember that any kernel for every quark 
sort traditionally carries the imaginary unit factor i 0, [|. The sign choice of A relies on this fact 
but, in principle, there is an alternative = A +L = —A L = —A +R for A > 0, if the kernels are 
redefined. An analysis shows that the denominator of Eq. (f|) is always positive and, therefore, the 
solution of this equation does exist at pretty large A. 

The quark matter state at finite chemical potential is defined by the saddle point of functional Eq. 
([!]) which we treat further maintaining the first nonvanishing contribution which is the figure-eight 
type diagram (see, for example, ||), 



1 — 2 (N c — 1) 



dp a(p) — 2 (3{p) 



7T 4 l-a(p) + /3(p) ' 

In our consideration this contribution to the generating functional (Z^ ~ e w ) would occur to be 

W = -np A InA + A (Y) , (Y) ~ I . (5) 
For the simplest situation of constant IL density the saddle point equation reads 

np A — A (Y) — . 

Comparing it with the gap equation Eq. @ one could notice the peculiar feature which is very 
practical to keep numerical calculations under control. The saddle point equation leads to the 
condition of gap independent of p. 

It was demonstrated in Refs. M , || that the quark back- reaction upon IL could be perturbatively 
estimated by studying the small variations of the IL parameters Sn and Sp around their equilibrium 
values n and p. Such variations could be incorporated by the IL theory if the deformable (crumpled) 
(anti-)instantons of size p, being the function of x and z, i.e. p — > p(x, z), are treated as the saturating 
configurations of the functional integral. Theoretically it is justified by developing the variational 
principle which makes use the Ritz method based on the multipole expansion of the deformation 
fields in the following form 

/x / \ dp 9 2 p 
p(x, z) = p(z, z) + — y u + . . . , y = x - z . (6) 

dz,j. p dZyflZv 

Then it is obvious the singular nature of the pseudoparticle solution generates the necessary cut-off 
at the scale of average instanton size. 

In principle, the instanton orientation in the colour space should be treated in similar manner. 
However, it does not lead to noticeable consequences while in the random instanton medium. Al- 
though such a treatment turns out very informative if one considers the polarization characteristics 
of the IL when the probe colour charge is present. 

In addition, the variations of zero modes in the interaction vertices of quark determinant at the 
transformation p —>■ p + 5p should be taken into account. Then as the output we have that for 
the long wave length excitations (for example, 7r-mesons) the deformation field describes colour less 
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scalar excitations of IL with the mass gap M of the order of several hundreds MeV, M 2 = ~ where 

v = - ~2 - , b = ~~~^~g ~ ) K is the kinetic coefficient being derived within the quasi-classical 
approach. N c and Nf are the numbers of quark colours and flavours, respectively. Our estimates 
give for this coefficient value of a few single instanton actions j3 = 8n 2 /g 2 , i.e. k ~ c (3 (with the 
factor c ~ 1.5 — 6 according to the ansatz taken for the saturating configurations) ||. 

Then, besides the diagrams with four legs (see, the term V in (p])) the extra diagrams with the 
scalar field attached (relatively speaking, the derivative in p of the vertex function in which the 
variation of the functions describing the zero mode 

dhi 

hi hi + -rr 1 $P , i = 1) • • ■ A , (7) 
op 

should be performed) are generated. Due to the diquark condensate presence handling the Lagrangian 
is substantially simplified if one restricts oneself with the leading contributions coming from the 
tadpole diagrams. Indeed, the leading contributions come from the term V and the term shown in 
Fig. 1 where two vertices are linked with the propagator of scalar field. 



Figure 1: The diagram of tadpole approach (see, the text). The solid (dashed) lines correspond to 
the fermion (scalar) field. 

Analyzing the modified Lagrangian it is easy to see that the form of the Gorkov equations would 
be the same if everywhere S is meant as an improved one £ + <5£ in which the correction <5£ is 
constructed with the modified functions Eq. (|7j). Moreover, the result for the Green function is 
formally the same if the modified functions E are meant. In particular, the form of gap equation is 
also retained because its kernel includes again the differences of the matrices which are constructed 
from the derivatives in p of matrices T, T + and their transposed forms. This difference is spanned 
by the unit vector as well. However, the following substitutions 

a(p) — * a(p) + Sa(p) , 

/3(p)-/?(p)+TO , 

where 

S a (p) = J — , 

are supposed to be performed and here J is the contribution of figure-eight type diagram with the 
scalar field propagator as an external leg 

J=^ 1_J_J 

iV^iVc-l) np^nAM 2 ' 
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Let us mention that effectively a dependence on the kinetic term k disappears (remember that 
|£) and its precise value is not operative in the approximation developed. The derivative of 



M 2 

the figure-eight type diagram which will shortly be necessary to proceed looks like 

1-0(P) * , n , a(p)-2 



8a(p) + 



(1 - a(p) + P(p)f 



And finally the modification of generating functional should be accommodated bringing Eq. (||) to 
the form 



W 



-4 

np 



nf )+X(Y). 



Then the new (if the variation of the IL density is absorbed) equation to calculate the saddle point 
reads 



-4 

np 



X (np )' In 
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np 
\X\ 



X(Y) = 



and the IL density is 
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T P ) r(z/ + i/2) 



1/2 



9 27 iV 9 • r 

= T .rf ?r is a measure of 



where the prime available means the derivative in A, the constant £ — 
interaction in the stochastic ensemble of PPs with an average size as 
pA = exp j- The derivative in A of the figure-eight type diagram looks like 




6(3 + a(3 + (3 2 r -4 + 7a + A(3 - ad 
— oa H 



a 



5/? 



and the derivative A' is defined by the following equation 



N C (N C 



2X 2 



2a; 

A 3 



dp 2(3 - 2f3 2 + 2a(3 - a 2 



7T 



(1-Q + /9)' 



Fig. 2 shows the results of calculating the parameter A (which is proportional to the free energy 
within the precision of one loop approximation) as the function of chemical potential at iV c = 3, Nf = 
2. The dashed lines expose a behaviour in the phase of nonzero values of diquark condensate (the 
left hand line (lower line) corresponds to the calculations if the quark interaction with IL is ignored) 
whereas the solid lines show the behaviours in the phase of broken chiral symmetry (upper line for the 
quark interaction with IL ignored) Pj. The saddle point parameter Ai of Ref. [0] and that exploited 



in the present paper are related as X\ 



X np 



The crossing points of the curves are fixing 



2(N C -1)N C 

the onset of colour superconducting phase. As mentioned above, actually, there is a transitional 
region of mixed phase where both chiral and diquark condensates develop nonzero magnitudes and 
the descent is not so steep. However, these details are inessential for further discussion. What is more 
interesting to be noticed comes from the crossing point of lower solid and dashed lines positioned on 
the plot. 

The first crossing point (along the /x-axis) which corresponds to the approach without the tadpole 
contributions included looks to be slightly larger in magnitude than that obtained earlier p c ~ 
300 MeV in Refs. |J, |J. However, it does not signal about shifting the critical point to larger values 
of p in the approach dealing with the initial Lagrangian. The IL parameters of various approaches are 
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Figure 2: The saddle point A as the function of chemical potential fx at N c = 3 and Nf = 2. The solid 
lines correspond to the phase of broken chiral symmetry (the lower one absorbs the IL perturbation). 
The dashed lines correspond to the colour superconductivity phase (the left hand one does not include 
the tadpole contribution whereas the right hand one does). 
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Figure 3: The quark density being defined by the crossing point of the upper solid and left hand 
dashed lines in Fig. 2. The solid line corresponds to the stable phase and the dashed one corresponds 
to the metastable one. 



sometimes slightly different and, in principle, could be optimized fitting, for example, Aqcd- Thus, 
the shift of fi c value pointed out is well within a precision of the IL theory. Regarding the second 
crossing point (when IL is perturbed by quarks) the conclusion could be more indicative. Here \x c 
becomes significantly larger. The corresponding quark matter densities as the functions of chemical 
potential for both approaches are depicted in Fig. 3 and Fig. 4. Apparently, on both plots the 
density values corresponding to the onset of colour superconductivity phase are appreciably larger 
than the density of normal nuclear matter. Perhaps, our result could be taken as a general indication 
on the considerable role of supplementary (insignificant on the instanton background) interactions of 
light quarks. They are able to change the estimate of critical density for colour diquark condensation 
drastically. 

Now we would like to comment on the crossing point of upper dashed and lower solid lines which 
is not shown in Fig. 2. It corresponds the unrealistic values of ft and in that region the IL is hardly 
applicable. One of the limits for IL approach coming from the very large values of chemical potential 
points out that with quark matter density increasing the average interquark distances may become 
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Figure 4: The quark density being denned by the crossing point of the lower solid and left hand dashed 
lines in Fig. 2. The solid line corresponds to the stable phase and the dashed one corresponds to the 
metastable one. 



very small and the 'Coulomb' (perturbative) field strengths would occur to be comparable with 
the (anti-)instanton ones. Then the (anti-)instanton superposition is not a proper configuration to 
saturate the functional integral. Thus, our conclusion resulting from the approach developed gives a 
message that the quark perturbation of IL leads the corresponding curves to move apart inherent in 
the phases of nonzero values of diquark and chiral condensates in Fig. 2. The 'chiral' curve becomes 
steeper but 'diquark' curve is displaced to the larger values of \x increasing its critical value. It looks 
like giving more reliability to our understanding of the perturbative fields role in the IL theory could 
provide us with more accurate appraisals for the phenomena considered. 
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Figure 5: The IL density for the approach at N c = 3 and Nf = 2. The solid line corresponds to the 
calculation in the phase of broken chiral symmetry. Both dashed lines correspond to the calculation 
in the phase of nonzero diquark condensate values. The lower dashed line gives the estimate of IL 
density without the tadpole contributions included. Regarding the upper dashed line see the text. 



Fig. 5 demonstrates the IL density as the function of /i when the quark back-reaction is incor- 
porated in the phase of broken chiral symmetry (solid curve) and in the phase of nonzero values of 
diquark condensate (dashed curve). It is interesting to notice that the curves manifest the different 
character of quark influence on IL. In the phase of broken chiral symmetry the density is almost con- 
stant in the suitable interval of /x whereas in the colour superconducting phase quickly disappears. 
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It is remarkable that for the former the corresponding analogue of the tadpole contribution (Jfijjj m 

Eq. (H) is strictly positive and, therefore, can lead to the increase of the IL density (gluon conden- 
sate) which is truly insignificant and demonstrates simply the approach sensitivity to the dynamical 
quark mass variation 0. In the latter case the sign of tadpole contribution is changing. As seen 
from Fig. 5 it leads to the drastic change of the IL density behaviour and the gluon condensate is 
getting weaker. Strictly speaking this effect has already been discussed in [1C]. 
The authors are partly supported by STCU #P015c, CERN-INTAS 2000-349, NATO 2000-PST.CLG 
977482 Grants. 
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